OPEN 3 ACCESS Freely available online 



•0-PLOS I ONE 



Patterns of Deep-Water Coral Diversity in the Caribbean (fi\ 
Basin and Adjacent Southern Waters: An Approach c^Ma* 
based on Records from the R/V Pillsbury Expeditions 

Ivan Hernandez-Avila 1 ' 2 * 

1 Departamento de Ciencias, Unidad de Cursos Basicos, Nucleo de Nueva Esparta, Universidad de Oriente, Margarita Island, Venezuela, 2IFREMER, Unite de recherche 
Etude des Ecosystemes Profonds, Laboratoire Environnement Profond, Plouzane, France 

Abstract 

The diversity of deep-water corals in the Caribbean Sea was studied using records from oceanographic expeditions 
performed by the R/V Pillsbury. Sampled stations were sorted according to broad depth ranges and ecoregions and were 
analyzed in terms of species accumulation curves, variance in the species composition and contributions to alpha, beta and 
gamma diversity. According to the analysis of species accumulation curves using the Chao2 estimator, more diversity occurs 
on the continental slope (200-2000 m depth) than on the upper continental shelf (60-200 m depth). In addition to the 
effect of depth sampling, differences in species composition related to depth ranges were detected. However, the 
differences between ecoregions are dependent on depth ranges, there were fewer differences among ecoregions on the 
continental slope than on the upper continental shelf. Indicator species for distinctness of ecoregions were, in general, 
Alcyonaria and Antipatharia for the upper continental shelf, but also the scleractinians Madracis myriabilis and Cladocora 
debilis. In the continental slope, the alcyonarian Placogorgia and the scleractinians Stephanocyathus and Fungiacyathus were 
important for the distinction of ecoregions. Beta diversity was the most important component of gamma diversity in the 
Caribbean Basin. The contribution of ecoregions to alpha, beta and gamma diversity differed with depth range. On the 
upper continental shelf, the Southern Caribbean ecoregion contributed substantially to all components of diversity. In 
contrast, the northern ecoregions contributed substantially to the diversity of the Continental Slope. Strategies for the 
conservation of deep-water coral diversity in the Caribbean Basin must consider the variation between ecoregions and 
depth ranges. 
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Introduction 

Deep-water corals represent a diverse group of cnidarians, 
including scleractinians, zoanthids, alcyonarians, antipatharians 
and hydrocorals (Stylasteridae), and are widely distributed in all 
oceans at depths below 50 m [1,2]. Despite its diversity, this group 
is ecologically important because it includes many habitat-building 
species (e.g., scleractinian formations and octocoral aggregations) 
and many species form associations with other species occurring 
on their surfaces [3], Coral formations represent one of the most 
important biotic habitats in the deep-water environment [4]. 
However, they are vulnerable to a range of anthropogenic threats 
that include trawl fisheries, oil exploitation, mining and ocean 
acidification [5-8]. 

The Caribbean Sea is a semi-enclosed basin of approximately 
2.75- 10 h km 2 bounded by Central and South America and the arc 
of the Greater and Lesser Antilles. The average depth is 
approximately 2400 m, with around 75% of the seafloor located 
at depths corresponding to abyssal waters and the continental rise 
(below 1 800 m) and the rest located on the continental shelf and 
continental slope. Hadal depths are uncommon and are restricted 
to the Cayman Trough and the Puerto Rico Trench. Sills and 



ridges separate the Caribbean into five basins: the Grenada, 
Venezuelan, Colombian and Yucatan Basins and the Cayman 
Trough [9], According to Spalding et al. [10], the Caribbean Sea 
includes five marine ecoregions: the Eastern Caribbean in the 
Grenada Basin; the Southern Caribbean in the Venezuelan Basin; 
the Southwestern Caribbean in the Colombian Basin; the Western 
Caribbean in the Yucatan Basin; and the Greater Antilles formed 
by Cuba, Hispaniola, Puerto Rico and the Cayman Islands. 

Topographic configuration of the Caribbean is very heteroge- 
neous in terms of depth, slope and elevations, such as volcanoes 
and seamounts. For example, more seamounts are found in the 
Greater Antilles and the Yucatan Peninsula and mud volcanoes 
have been found in the Southern Caribbean [9,1 1,20]. In general, 
deep-water hydrography of the Caribbean consists in stratified 
masses of water in the upper 2000 m and unstratified water below 
this depth [12]. Below the mixed layer are found Subtropical 
Underwater (SUW 110-200 m), composed mainly of Southern 
SUW and Eastern Tropical Water, both entering through passages 
in the Lesser Antilles, and Northern SUW that inflows through the 
Windward Passage and the Anegada-Jungfern complex. Central 
Water (CW), composed mainly of eastern South Atlantic CW and 
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secondarily of North Atlantic CW, follows a pattern of inflow 
similar to the SUW. Antarctic Intermediate Water (AAIW) is the 
main water mass between 800-1 100 m and inflows mainly via the 
Grenada Passage. Deep inflow of North Atlantic Deep Water and 
eventual outflow occur between 1200-2000 m through the 
Windward Passage and the Anegada-Jungfern complex [13,14]. 
Below 2000 m, elevated sills limit the interchange between the 
Caribbean and the Adantic, although episodic ventilation could 
occur [14]. The Yucatan Peninsula serves as a gateway for most of 
the outflow water from the Caribbean [15]. Additionally, the flux 
of particulate organic carbon (POC flux) shows spatial differences 
in the Caribbean. Models [16,17] estimate more POC flux to the 
deep benthos in the Southern Caribbean due to the influence of 
coastal upwelling and riverine inflow. Moreover, the Lesser 
Antilles could be affected by the Orinoco-Amazonas inflow 
[18,19]. 

Deep-water corals in the Caribbean have been studied 
frequently since the 19 th century (see historical review in [20]). 
Although most of these studies have focused on taxonomy and 
systematics, several studies have described the large-scale distri- 
bution of deep-water corals in the Adantic [21-23,34] or at the 
meso- or local scale in certain regions of the Caribbean [24—26] . 
These studies highlight the Caribbean as a globally important 
center of diversity for deep-water corals [1] . Most of the margins of 
the Caribbean Basin offer suitable habitats for deep-water corals 
[27-29], and many habitat-building species are widespread, 
forming deep-water reefs and similar habitats [20,24]. However, 
the general patterns of deep-water coral diversity in the Caribbean 
remain elusive, as is the case for most of the deep-water benthos 
due to fragmented and (apparently) insufficient sampling efforts 
[9]. 

Several models have been proposed to describe general patterns 
of deep-water diversity relative to depth and spatial distributions. 
In an analysis of marine diversity according to depth, Gray et al. 
[30] postulated that shallow and deep habitats support similar 
numbers of species in similar areas and that the overall diversity 
relative to the number of habitats and number of species per 
habitat must be higher in the shallow-water benthos than in the 
deep-water benthos. Other studies support the notion that there is 
higher species diversity in assemblages of the deep-water benthos, 
in relation to assemblages of the continental shelf, for particular 
taxonomic groups and for complete community assemblages 
[7,31,32]. For deep-water corals, Cairns [1] has shown that most 
species inhabit the depths of the continental slope, although the 
ratio of species inhabiting the shelf (depths of 50-200 m) to species 
inhabiting the slope is relatively high. 

Many deep-water coral species are widely distributed, and 
models of habitat suitability usually include broad connections 
between regions [28,29] . These patterns seem to indicate that the 
composition of deep-water corals could be homogeneous at a 
regional scale. For the Caribbean Sea, Miloslavich et al. [9] failed 
to detect regional differences in shallow-water faunal composition 
for five major taxonomic groups, including corals. This result 
suggests that the classification of marine ecoregions in the 
Caribbean did not apply to shallow-water benthic fauna. 
However, according to Levin et al. [32], deep-water species 
diversity shows geographic variation at a scale of 100-1000 km, 
the same scale as that of ecoregions in the Caribbean, and Briggs et 
al. [33] report regional differences in species composition in the 
abyssal depths of the Caribbean. Furthermore, the diversity 
contours of azooxanthellate Scleractinia presented in Cairns [1] 
and the cluster analysis in Cairns & Chapman [23] suggests that 
species are not homogeneously distributed in the Caribbean Basin. 
In contrast, Dawson [34] re-analyzed previous data [21-23] and 



did not find evidence of subdivision in regions within the 
Caribbean based on the distribution of species, excluding 
surrounding regions, although an effect of depth in species 
composition was detected. Therefore, there is an on-going debate 
on whether the spatial distribution of deep-water corals in the 
Caribbean region is heterogeneous [23] or homogeneous [34], 
which is very important to resolve as it will determine the most 
appropriate conservation strategy for this ecosystem. 

The identification of the patterns of spatial distribution of the 
deep-water benthos on a large scale is limited by the costs and 
logistic challenges involved in studying these systems. Also, 
analyses of databases are usually difficult due to differences in 
the time of sample collection, methods, identification and 
coverage. Previous approaches related to distribution of deep- 
water corals in the area [23,34] were based on ordination methods 
of data pooled into subunits, which offers an excellent framework 
but overlooks the within-unit variation of assemblages in order to 
detect differences among units. However, a sampling program 
specific to the Caribbean could provide the opportunity to conduct 
an analysis of general patterns. 

An extensive sampling program in the Caribbean and adjacent 
waters was conducted by the R/V Pillsbury of the Rosenstiel 
School of Marine and Atmospheric Science (RSMAS), University 
of Miami. These samples were collected at approximately 800 
stations in shallow- and deep-water habitats over a relatively short 
period (approximately 5 years). These expeditions produced 
records of corals below a depth of 60 m for 204 stations in the 
Caribbean Basin and Guiana (Fig. 1), and were assembled 
together in a database that provides a unique opportunity to 
overcome the problems stated above regarding the identification 
of spatial distribution patterns in the Caribbean. Although the 
dataset from the Pillsbury expedition does not comprise all species 
and records of deep-water corals in the Caribbean, the number of 
species and records and the distribution of stations suggest that the 
dataset could furnish a representative sample for determining 
general patterns of diversity. The aim of this paper is to determine 
the patterns of spatial distribution of deep-water coral diversity in 
the Caribbean in terms of wide depth ranges and ecoregions, 
including the adjacent Guianian ecoregion, based on the records 
of the Pillsbury expeditions. 
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Figure 1. Stations of deep-water coral records analyzed in the 
present study. EC: Eastern Caribbean, Gui: Guianian, SC: Southern 
Caribbean, SWC: Southwestern Caribbean, WC: Western Caribbean. 
Regions adjusted from Spalding et al. [10]. 
doi:10.1371/journal.pone.0092834.g001 
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Materials and Methods 

The database contains records of deep-water corals collected by 
the R/V Pillsbury from 1966 through 1971 in the Caribbean and 
the adjacent southern waters (Guianian region). The complete 
depth range covered between approximately 60 and 5400 m, 
representing 204 stations bordering almost the entire Caribbean 
basin (Fig. 1), except for part of the Central American region 
where the shelf is shallow. Samples were collected with 10-ft or 41- 
ft otter trawls. Specific information for each sample is found in 
Voss [35-38], Staiger [39-41], and Staiger & Voss [42]. 
Specimens were identified by major coral specialists (FM Bayer, 
SD Cairns, DM Opresko and collaborators) and deposited in the 
U.S. National Museum of Natural History (Smithsonian Institu- 
tion-NMNH) and the Marine Invertebrate Museum (MIM) of 
RSMAS, University of Miami. Station and depth data related to 
the specimens deposited in the NMNH was obtained from the 
collection database (http:/ /invertebrates. si.edu/collections.htm), 
and for the species held in MIM-RSMAS by examination of 
records. Taxonomic records were revised according to Williams 
and Cairns (http://researcharchive.calacademy.org/research/ 
izg/OCTOCLASS.htm), the World Marine Species Database 
(www.marinespecies.org/) and taxonomic treatises [21,22]. Any 
names not found in these sources were excluded from the analysis. 
The database was sorted by presence-absence of species by station, 
avoiding duplicate records among data sources. 

The sampling stations were classified according to Spalding et al. 
[10] based on the following ecoregions: Greater Antilles (50 
stations,st), Eastern Caribbean (35 st), Southern Caribbean (45 st), 
Southwestern Caribbean (43 st), Western Caribbean (14 st) and 
Guianian (17 st). These ecoregions comprise the border of the 
Caribbean Basin and its adjacent southern ecoregion, which were 
divided into sub-basins and defined by consensus of biogeographic 
studies and expert opinions [10]. This classification is similar to the 
regions proposed by Cairns & Chapman [23] for deep-water 
Scleractinia deduced by cluster analysis (except by division of the 
"insular" cluster [23]). Nevertheless, the northern coasts of Cuba, 
Hispaniola and Puerto Rico, assigned by Spalding et al. [10] to the 
Bahamian ecoregion, are considered part of the Greater Antilles 
for the purpose of this study. Depth ranges were separated 
according to McClain & Hardy [43] except for the inclusion of the 
upper continental shelf, ranging from 60 to 200 m (103 St.). Other 
depth ranges considered include the continental slope (200- 
2000 m- 93 St.), the continental rise (2000-4000 m- 5 st.) and the 
abyssal depths (4000-6000 m- 3 St.). The boundaries of the upper 
continental shelf and the continental slope in this classification 
coincide with the common criteria for defining shallow and deep 
corals habitatsfl, 2 1-23,34,44] and the low occurrence of Scler- 
actinia below 2000 m depth [1]. The other depth ranges were 
selected based on commonly accepted marine depth ranges [10]. 

To illustrate the patterns of spatial distribution of species relative 
to depth, the distribution range of each species was estimated 
according to the minimum and maximum depths recorded in this 
study. Moreover, for the upper continental shelf and the 
continental slope, species accumulation curves were plotted using 
two approaches: 1) the accumulation of species by sampling effort 
and 2) the Chao2 estimator [45,46]. The Chao2 estimator [45] 
was selected for this analysis because it is a non-parametric 
approach that considers the frequency of each species in the 
samples. This method shows high precision for estimating 
simulated communities [47] and has been used for the analysis 
of the richness of marine fauna [47-49]. Although the Chao2 
estimator is more conservative than some parametric methods 
[50], it allows the total species richness to be compared between 



depth ranges without parametric assumptions [51]. To avoid the 
bias resulting from the selection of only stations with positive 
records, the estimates of coral richness were derived from 
calculations that also included stations with negative coral records 
from the expedition but those that included other megafaunal 
invertebrates (e.g., echinoderms, gastropods, decapod crusta- 
ceans). Differences in species accumulation curves between depth 
ranges were studied by graphical comparison of their respective 
95% confidence intervals [52], and by t-student test following 
Guerra-Castro [53]. 

The null hypothesis of no differences among depth ranges and 
regions was tested using a permutational multivariate analysis of 
variance (PERMANOVA, [54]) on a matrix of Sorensen's 
similarity index for each pair of samples. The depth of each 
sample was included as a covariable in the analysis. Variation 
among samples within the same ecoregion and depth range was 
used to generate the residuals for the model. The distances 
between samples according to their similarity were represented in 
a non-metric multidimensional scaling plot (MDS). The factor 
Depth Range had only two levels: upper continental shelf and 
continental slope. Records for the continental rise and abyssal 
ranges were scarce and did not allow comparisons. Additionally, 
permutational analysis of multivariate dispersions (PERMDISP) 
[55] and pairwise PERMANOVA tests between regions were 
performed for each depth range. Frequent (i.e., indicator) species 
contributing to differences between ecoregions and depth ranges 
were identified using an information statistic (/-test, 2AI) [56]. 

The relative contribution of the regions to the overall diversity 
was estimated based on the additive partitioning of diversity 
according to Lu et al. [57]. For this purpose, species lists were 
pooled by ecoregion, grouped according to the pairwise test for the 
upper continental shelf and for the continental slope (except for 
distant ecoregions), and depth range. Gamma diversity, repre- 
senting the total richness recorded in the present study, was 
estimated as an additive partitioning of alpha and beta diversity: 
y = a + /?.The overall a diversity, the richness for each ecoregion 

1 s 

and depth combination, was expressed as a s t= «,, where n is 

i 

the number of ecoregion-depth pairs, and ni the number of 
ecoregions with records of the i th species. The partitioning of a for 
the k th ecoregion-depth combination was calculated as 
* j 

0L s (k) = > -, where s k represents each record for the ecoregion- 
I 

depth combination. Similarly, beta diversity was calculated with 

Sk Sk / \ 

n T^n — ili .. . y-^(n — Hi) 

pj= > , partitioned using p (K)= > , and gam- 

■ L — < n z — 1 n-n, 

i i 1 

ma diversity with y s t = a s t + f} s t, partitioned using 

Sk J 

y s (k) = — = a. s (k) + fi s (k). The components of diversity were 
i n ' 

estimated for each ecoregion to evaluate their contribution to 
gamma diversity. In this case, this approach was not used to 
estimate the overall diversity, but to estimate the relative 
importance of alpha and beta diversity and the relative contribu- 
tion of each ecoregion-depth combination. An estimate of 
between ecoregion-depth differentiation was calculated as 

D s t= — . A value of D„>0.5 means that most diversity is 

distributed among units (ecoregion-depth combinations) [57]. 
After estimating the components of diversity partitioned by 
ecoregion, they were expressed as relative contribution to 
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In certain cases, taxa were identified only to the generic level. In 
these cases, the taxon was assumed to differ from other species 
identified in the same genus. However, if identification at the 
generic level corresponded to a species previously recorded in the 
database, the partitioning of records could affect the data analysis. 
To estimate if this effect influenced the results, all data were 
aggregated to the generic level, and certain data analyses, 
including taxa accumulation according to sampling effort, the 
PERMANOVA test, the pairwise test and the MDS plot, were 
repeated and compared with the analysis based on identification to 
the species level. In addition, a RELATE analysis [55] was 
performed to compare the matrices of resemblance between 
samples based on identification to the species level and to the 
generic level. The results of these analyses are shown in the Tables 
SI in File SI and Table S2, Fig. S1-S2 in File S2. 

Results 

A total of 2 1 8 species were recorded, belonging to almost all 
major groups of deep-water corals (105 Alcyonacea, 38 Anti- 
patharia, 73 Scleractinia, 2 Anthoathecatae). The dataset includes 
935 records from 204 stations and both colonial and solitary coral 
species. Alcyonarians, antipatharians and azooxanthellate scler- 
actinians comprised most of the species, but several zooxanthellate 
and aposymbiotic species were recorded, especially from the upper 
continental shelf (Table S3 in File S2). Although some taxa have 
few records and species compared with previous records for the 
Caribbean Basin, such as Stylasteridae [58,59], in general the 
dataset is a representative sample of the deep-water coral diversity 
of the area. 

Most of the species were found at depth intervals that included 
the continental slope. The next most common group consisted of 
the species found on the upper continental shelf (Fig. 2). In all, 
58.2% of the species whose depth intervals included the upper 
continental shelf were also present below this depth range. A total 
of 62.9% of the species whose depth intervals included the 
continental slope were restricted to this depth range. More species 
(151 species) were found on the continental slope than on the 
upper continental shelf (123 species). The continental rise and 



abyssal ranges harbored few species (7 and 4 species, respectively), 
but the low sampling effort in these ranges precluded any 
comparative approach. 

The species accumulation curve of the upper continental shelf 
estimated a richness of 123±4.69 species (mean ± sd) (n= 138), 
and 151 ±8.44 species for the continental slope (n=138) 
(t 9 99 = 91.65, p<0.0001). The Chao2 estimator for species on 
the upper continental shelf was 148. 02± 10.98 species, slightly 
more than one-half of the species predicted for the continental 
slope (279.18±38.81 species) (t 999 = 102.09, p<0.0001). The 
shape of the accumulation curve for the continental slope 
suggested that the species estimate could increase with more 
sampling. For either estimator, the 95% confidence intervals of the 
curves of upper continental shelf and continental slope did not 
overlap at the end of the curves (Fig. 3). 

The PERMANOVA test based on Sorensen's similarity index 
between samples indicated significant differences in coral assem- 
blages related to the depth of the sample, depth ranges and 
ecoregions (Table 1). A significant interaction was found between 
depth ranges and ecoregions, which suggest that ecoregional 
patterns depended on depth ranges. Differences between the 
composition of the upper continental shelf and continental slope 
species were detected in all ecoregions (pairwise test, p<0.001 in 
each case) except for the Western Caribbean due to a lack of 
samples from the upper continental shelf. A non-metric MDS 
based on the resemblance among all stations produced a 
degenerate solution (an accumulation of rankings of similarities 
in extreme positions of the plot). This outcome was associated with 
the occurrence of few stations that did not share their species 
composition with any other station due to the extreme rarity of the 
species (e.g., stations with only one species that did not occur in the 
other samples had a similarity of 0). These cases represented 
approximately 6.4% of the studied stations and 1.6% of the species 
records. A non-metric MDS plot excluding these stations showed 
differences in composition between the upper continental shelf and 
the continental slope (Fig. 4A). /-tests between depth ranges 
showed that indicator species of the upper continental shelf were 
mainly alcyonarians and antipatharians, but for the continental 
slope, most indicator species belonged to order Scleractinia 
(Table 2). Although the continental slope had higher species 
richness, fewer indicator species were found there than on the 
upper continental shelf, suggesting that differences were sustained 
in part by infrequent (rare) species. The continental slope showed 
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Figure 2. Depth intervals of species according to sample. Dashed lines represent limits of depth ranges. 
doi:1 0.1 371 /journal.pone.0092834.g002 
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Table 1. PERMANOVA for similarity in species composition of 
deep-water corals (Sorensen). Bold denote significant 
differences. 
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more rare species (i.e., single record) than the upper continental 
shelf (% 2 — 25.37, p<0.001, df= 1) although sampling efforts were 
similar. 

According to the a posteriori test, the similarity among ecoregions 
changed according to depth range (Table 3, Fig. 4B). For the 
upper continental shelf, differences were not related to changes in 
multivariate dispersions between ecoregions (PERMDISP, 
F=1.10, p = 0.90, df d = 95), which support the detection of 
species turnover between ecoregions. Differences in assemblage 
composition were found in all ecoregions except for the Greater 
Antilles-Eastern Caribbean group, and the Guianian ecoregion, 
which showed p-values that did not allow the assumption of 
differences or similarities. Ecoregions on the continental slope 
showed differences in multivariate dispersions (PERMDISP 
F = 4.138, p = 0.431, df d = 86) due to low dispersion of the 
Guianian ecoregion. Excluding the former ecoregion, the multi- 
variate dispersions were not different (PERMDISP F = 2.439, 
p = 0.125, dfd = 82). The Greater Antilles was different from the 
Eastern Caribbean and Western Caribbean. These latter two 
ecoregions also were different from the Southern Caribbean- 
Southwestern Caribbean group (Fig. 5). However, no differences 
were found between the Greater Antilles and Southern ecoregions, 
that shown an overlap of assemblages between separate regions. 

Differences among ecoregions were explained in part by 
changes in frequency of occurrence of coral species as indicated 
by the /-test (Table 4). For the upper continental shelf, the Greater 
Antilles-Eastern Caribbean group differs from others ecoregions 
because it has higher frequencies of Madracis myriaster and Micella 
guadalupensis. The Southern Caribbean ecoregion showed high 
frequency of the alcyonarians Bellonella rubistella, Aphanipathes pedata 
and Thesea sp. and the scleractinians Balanophyllia dineta and 
Cladocora debilis. The Southwestern Caribbean ecoregion showed a 
higher frequency of Diodogorgia sp. than the Greater Antilles. The 
Guianian did not show indicator species, which is consistent with 
the pairwise analysis. For the continental slope, the scleractinians 
were more important to determine the differences between 
ecoregions. In the Greater Antilles ecoregion, Stephanocyathus 
diadema was important, and was also frequent in the Southern 
Caribbean-Southwestern Caribbean along with Stephanocyathus 
paliferus. Higher frequency of Fungiacyathus symetricus was important 
for differences in the Eastern Caribbean ecoregion. The Western 
Caribbean ecoregion had no species that showed more frequencies 
than in other ecoregions. 

The differences in species composition among ecoregions 
suggest that the diversity of corals in the Caribbean Sea is a 
product of a high level of species turnover between ecoregions. 
This hypothesis is supported by the analysis of the additive 
partitioning of overall diversity (gamma diversity) into local 
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Figure 3. Accumulative species curves (S obs and Chao 2) for 
the upper continental shelf (red) and the continental slope 
(blue), dashed lines represent 95% confidence intervals. 

doi:10.1371/joumal.pone.0092834.g003 

diversity (alpha diversity) and turnover (beta diversity) [57]. Beta 
diversity was the most important component of gamma diversity 
(D st = 0.752). Moreover, ecoregions varied in their contribution to 
gamma and beta diversity. For the upper continental shelf, the 
Southern Caribbean showed the greatest contribution to overall 
alpha, beta and gamma diversity, followed by the Greater Antilles- 
Eastern Caribbean group. On the continental slope, the Greater 
Antilles, Western Caribbean and Eastern Caribbean ecoregions 
had important contributions to all components of diversity, and 
together they represented 41.1% of the gamma diversity (Fig. 6). 
The contribution of the Guianian ecoregion and the Southwest- 
ern-Southern Caribbean group were lower than other ecoregions 
in this depth range. 

Discussion 

Based on analyses of species accumulation curves, the present 
study confirms previous suggestions, based on species depth 
ranges, that coral diversity is higher on the continental slope in the 
Caribbean Sea [34]. Furthermore, the present study provided a 
robust test of these differences given the sampling coverage of the 
Pillsbury expeditions. An increase in faunal richness in the deep 
waters of the continental slope or rise relative to continental shelf 
habitats has been found by numerous studies [7,32,34,60-65], 
although these differences depend on the basins, taxa and types of 
communities studied. The critical views and evidence furnished by 
Gray et al. [30] show exceptions to this trend and emphasize the 
effect of the area sampled on some estimates. For deep-water 
corals, niche-specific boundaries could contribute to the observed 
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Figure 4. A. nm-MDS of the similarity between stations, 
estimated using Sorensen's similarity index, according to 
depth range. B. nm-MDS based on distance between centroids of 
ecoregions. Red, upper continental shelf; blue, continental slope. GA: 
Greater Antilles, EC: Eastern Caribbean, WC: Western Caribbean, Gui: 
Guianian, SC: Southern Caribbean, SWC: Southwestern Caribbean. 4A, 
12 outliers were excluded. 4B, WC from the upper continental shelf 
excluded due lack of samples. 
doi:1 0.1 371 /journal.pone.0092834.g004 

diversity pattern. According to Cairns [1], 63.4% of azooxanthel- 
late Scleractinia live between 200-1000 m, 14.7% live between 
1000-2000 m and 52.9% live between 50-200 m (these percent- 
ages sum to a value greater than 100% because the original values 
included species that inhabit more than one of the specified depth 
ranges). 

Although the diversity of the continental rise and abyssal depths 
was not estimated due to insufficient sampling, few coral species 
were present below 2000 m (Fig. 2). In these depth ranges, data 
records indicate that the species present were antipatharians 
(Antipathes umbratica, Abyssopathes lyra, Bathypathes seculata, Bathypathes 
sp. and Cirrhipathes sp.), scleractinians (Deltocyathus italicus, Fungia- 
cyathus cripus, F. marenzelleri) and one alcyonarian (Thesea sp.). 
Despite, Briggs et al. [33] did not report any deep-water coral 
species in samples from the Venezuelan Basin at depths between 
3411 and 5062 m [33], Pillsbmy data records indicated the 
presence of Thesea sp., F. cripus and A. lyra at about 4200, 4500 
and 5400 m, respectively. According to Guinotte et al. [6] the 
aragonite saturation horizon (ASH) in the waters of the 
northwestern Adantic is located between 2000-3000 m, limiting 
(although not prohibiting) the development of carbonate structures 
below this zone. Likewise, globally, only 4.7% of the deep-water 



scleractinians are found at depths below 2000 m [1] and global 
models suggest that the most suitable habitats for framework- 
forming corals and octocorals worldwide are the continental shelf 
and the continental slope [28,29]. Based on the current data, a 
ummodal curve can be used to describe the diversity of deep-water 
corals, with high values of richness on the continental slope and 
lower values in shallower waters and at bathyal depths, which is 
consistent with the pattern proposed by Dawson [34] for 
azooxanthellate Scleractinia in the Caribbean and surrounding 
regions based on depth ranges of the species. 

This study detected a turnover of coral species composition 
associated with large depth ranges, in addition to the effect of 
depth. Zonation by depth is a widespread pattern in the deep- 
water biota [66] . The overall diversity (gamma diversity) of deep- 
water corals depends on the number of species at each depth 
interval and the rate of species replacement between depth 
intervals. Changes in composition along depth ranges could be 
related to differences in habitat characteristics rather than depth 
between the upper continental shelf and the continental slope. 
Despite its depth, the continental slope shows high topographic 
and hydrographic heterogeneity [32,67], variation in slope 
inclination and seafioor, habitats (ridges, seamounts and table- 
mounts [9,20,27]), as well as in tectonic, terrestrial and oceanic 
influences that shape physical configurations, which have different 
hydrographic conditions such as water mass, temperature and 
carbon flux. The Caribbean continental slope is bathed by at least 
four water masses at different depths and points of entrance in the 
basin [12,13]. Furthermore, because both topographic and 
hydrographic conditions have an influence in the distribution of 



Table 2. Indicator species for each depth range as 
determined using the /-test. 



Species 


UCS (101) 


CS (92) 


2AI 


Acanella sp. 


0 


26 


18.68 


Diodogorgia sp. 


27 


0 


16.92 


Ellisella sp. 


31 


1 


16.41 


Stephanocyathus diadema 


0 


19 


13.22 


Fungiacyathus symmetricus 


0 


13 


8.82 


Bellonella rubistella 


14 


0 


8.31 


Antipathes lenta 


24 


3 


8.27 


Deltocyathus eccentricus 


0 


11 


7.4 


Deltocyathus italicus 


0 


10 


6.7 


Leptogorgia sp. 


11 


0 


6.45 


Swiftia excerta 


11 


0 


6.45 


Tanacetipathes barbadensis 


10 


0 


5.84 


Balanophyllia dineta 


10 


0 


5.84 


Nicella goreaui 


9 


0 


5.24 


Nicella guadalupensis 


9 


0 


5.24 


Scleracis sp. 


15 


2 


4.74 


Telesto sp. 


12 


1 


4.59 


Cladocora debilis 


7 


0 


4.04 


Stephanocyathus coronatus 


0 


6 


3.95 


Stephanocyathus paliferus 


0 


6 


3.95 



Occurrence in samples of upper continental shelf (UCS) samples and in 
continental slope (CS) samples are given. The total sample size is given in 
parentheses. Only species with 2AI>3.84 are shown. 
doi:1 0.1 371 /joumal.pone.0092834.t002 
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Table 3. Type I error probabilities for pairwise f-tests on the similarity in deep-water coral species composition between 
ecoregions (columns vs. rows). 





Upper continental shelf 




GA 


EC 


Gui 


SC 


SWC 


EC 


0.1322 










Gui 


0.0009 


0.0330 








SC 


0.0001 


0.0001 


0.0379 






SWC 


0.0001 


0.0019 


0.0215 


0.0001 




WC 


0.3366 


0.1142 


0.2547 


0.1532 


0.2714 


Continental slope 




GA 


EC 


WC 


SWC 


SC 


EC 


0.0072 










WC 


0.0143* 


0.1627 








SWC 


0.0429 


0.0001 


0.0015 






SC 


0.1804 


0.0032 


0.0016 


0.1302 




Gui + 


0.0757 


0.001 1 


0.0063 


0.1780 


0.6011 



GA: Greater Antilles, EC: Eastern Caribbean, Gui: Guianian, SC: Southern Caribbean, SWC: Southwestern Caribbean, WC: Western Caribbean. Significant differences are 
shown in bold. Gray denote low sampling effort. 
"Marginal p-value. 
+PERDISM p<0.05. 

doi:1 0.1 371 /journal.pone.0092834.t003 



deep-water corals [28,29,67,68], heterogeneity in these conditions 
could enhance fauna diversity [69,70]. 

The differences detected between ecoregions in species compo- 
sition on both the upper continental shelf and the continental slope 
support the hypothesis that the deep-water fauna of the Caribbean 
Basin and nearby waters reflect a heterogeneous species compo- 
sition at a regional scale. This finding contrasts with previous work 
[34] that suggests the absence of differences in composition of 
deep-water corals between ecoregions in the Caribbean basin, for 
both depth ranges. Nevertheless, similarity between nearby 
ecoregions, forming larger groups (i.e., Greater Antilles-+Eastern 
Caribbean on the continental shelf and the Southwestern 
Caribbean-Southern Caribbean group), and between distant 
ecoregions (i.e., Greater Antilles and Southern ecoregions on the 
continental slope) suggests a wide connectivity at both depth 



ranges. Gene flow and connectivity of deep-water corals has been 
associated with water mass stratification [11,12]. Henry [2] 
postulated that widespread circulation of the Antarctic Interme- 
diate Water in the Caribbean could enhance larval flow of deep- 
water corals through the basin. 

The ecoregions obtained by the present study differs from those 
proposed by Cairns & Chapman [23] for deep-water Scleractinia; 
however both studies show heterogeneity in the Caribbean Basin, 
and affinities between ecoregions. The discrepancies related to 
ecoregions between the previous works [23,34] and the present 
study mainly involve detection of within-ecoregion variation (in 
contrast with Dawson [34]), and the separation of the northern 
regions, sorted previously [23] as "insular" cluster, which could be 
explained by the increase of taxa coverture (by inclusion of 
Alcyonaria and Antipatharia), data source (use of samples from 




Figure 5. Interpretation of pairwise analysis between ecoregions by depth ranges. (A) Upper continental shelf, (B) continental slope. 
Affinities with the Guianian ecoregion were not resolved (see Results) and affinities between the Greater Antilles and Southern ecoregions for the 
continental slope are not shown. 
doi:1 0.1 371 /journal.pone.0092834.g005 
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Table 4. Indicator species in each ecoregion detected using the /-test. 







Upper continental shelf 










GA+EC 


SC 


SWC 


Gui 


GA+EC 




Bellonella rubistella 


Diodogorgia sp. 




Balanophyllia dineto 


Leptogorgia sp. 


Aphanipathes pedata 


Cladocora debilis 


SC 


Madracis myriaster 








swc 


Madracis myriaster 


Thesea sp. 








Nicella guadalupensis 


Balanophyllia dineta 






Aphanipathes pedata 


Cladocora debilis 


Bellonella rubistella 


Gui 


Madracis myriaster 


Bellonella rubistella 






Thesea sp. 


Continental slope 




GA 


EC 


swc+sc 


WC 


GA 




Fungiacyathus symmetricus 


Stephanocyathus paliferus 




EC 


Stephanocyathus diadema 




Stephanocyathus diadema 




WSC+SC 




Fungiacyathus symmetricus 






WC 






Stephanocyathus diadema 





Ecoregions were grouped according to pairwise analysis. GA: Greater Antilles, EC: Eastern Caribbean: WS: Western Caribbean, SWC: Southwestern Caribbean, SC: 
Southern Caribbean, Gui: Guianian. 

Pairwise comparisons are to be read as columns vs. rows: species are more frequent in the ecoregion (or group) heading a column than in the ecoregion of the 
intersecting row. Only species with 2AI>3.84 are shown and those with values of 2AI>6.63 are shown in bold. 
doi:1 0.1 371 /joumal.pone.0092834.t004 



one expedition survey) and statistical approach (analysis of 
variance instead of hierarchical clustering). Although a cluster 
analysis of checklist could be preferable by the selection of a 
"complete" list of species, using records from different surveys 
could introduce selection bias due to the area and ecoregion(s) 
covered in each survey and due to pooling data without estimating 
the variation between samples within the same unit (i.e., 
ecoregions). Moreover, data manipulation, such as excluding 
endemic species [34], could overlook an important component of 
diversity. Although endemism could be an artifact of sampling 
coverage and collection methods [71], excluding them may 
actually contribute to the homogeneity previously reported for 
the Caribbean basin [34]. 

Despite differences in sampling design or checklist studies, the 
analysis of indicator species (/-test, Table 4) was generally 
consistent with other studies. On the upper continental shelf, 
Madracis myriabilis has widespread distribution in the Caribbean 
but is more frequent in the Antillean arc [21], although there are 
exceptions [21,22,25]. Micella guadalupensis have a similar pattern of 
distribution; there are few records of this species in offshore reefs 
from Colombia [24,44] and Venezuela [72], compared to its 
distribution along the Antilles [72]. There are species that are 
frequent in the Southern Caribbean but not in the Southwestern 
Caribbean, such as Aphanipathes pedata, which has not been 
reported in Colombian Shelf waters to date [44]. Similarly, 
Balanophyllia dineta and Bellonella rubistella have been reported with a 
frequency lower than 5% for samples from Colombia [44], 
compared to 40-43% from the Southern Caribbean in the present 
study. The occurrence of Cladocora debilis in the Southern 



Caribbean in contrast with the Greater Antilles is consistent with 
the distribution pattern of this species [22], but differences in 
frequency between the Southern and Southwestern Caribbean 
have not been reflected in recent samples [25,44]. Other taxa such 
as Leptogorgia sp. and Thesea sp. cannot be compared because they 
contain various species with different patterns of distribution. 
Regarding the frequent species in the Southwestern Caribbean, 
Diodogorgia sp., could be related with the only species of the genus 
reported, Diodogorgia nodulifera, which is a widespread species on the 
Colombian shelf [44]. 

Although the continental slope had higher diversity, only a few 
species showed contrasting frequencies, in part because the 
assemblages showed less variation between ecoregions at this 
depth (Table 3), and the variation could be enhanced by 
infrequent species. The frequency of Stephanocyathus species (S. 
paliferus and S. diadema) in the Southwestern Caribbean-Southern 
Caribbean group is supported by recent samples in Colombian 
deep waters [44] , but both species show widespread distribution in 
the Caribbean [21]. However, Fungicyathus symetricus was found 
only in low frequencies in the Southern and Southwestern 
Caribbean ecoregions, being recorded only in the deep waters of 
offshore Colombian islands [44], compared to numerous records 
in the Eastern Caribbean [21]. 

The substantial contribution of beta diversity to overall diversity 
indicates that the high diversity of deep-water corals is a product of 
high turnover in assemblages between ecoregions and depth 
ranges. A potential explanation is that the observed turnover is 
associated with the interaction between environmental variation 
(topographic and hydrographic) and niche-specific conditions. 
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Figure 6. Relative contribution of regions to the diversity of 
deep-water corals for the upper continental shelf (red) and 
continental slope (blue) in the Caribbean Basin. EC: Eastern 
Caribbean, Gui: Guianian, SC: Southern Caribbean, SWC: Southwestern 
Caribbean, WC: Western Caribbean. 
doi:1 0.1 371 /journal.pone.0092834.g006 

Differences in deep-water fauna at global and regional scales have 
been related to environmental factors, including productivity, 
sediment heterogeneity, oxygen availability, hydrodynamics, slope 
and disturbance [28,29,32,73-76]. In addition to the topographic 
variation, the Caribbean shows differences in POC flux [16,18] 
and other sources of carbon and energy in the Caribbean can 
contribute direcdy to the benthos. These sources include 
hydrothermal vents, hydrocarbon seeps and mud volcanoes, all 
of which could affect faunal assemblages [77-79]. An approach to 
the hypothesis of environmental factors as causes of regional 
differences in coral composition requires correlative patterns of 
assemblages with environmental models, for both water mass 
characteristics and topographical attributes. Future work should 
attempt to determine if the patterns detected in corals are found in 
other megafaunal invertebrates (crustaceans, mollusks, echino- 
derms, etc.) from deep waters. 

The importance of inter-regional beta diversity and the relative 
contributions of different ecoregions to gamma diversity highlight 
the need for joint efforts by different countries to facilitate the 
management and conservation of deep-water diversity. In general, 



the Caribbean Basin and nearby regions face potential direct 
anthropogenic impacts on the deep-water benthos produced by 
bottom trawl fisheries and oil exploitation [8,80,81]. Moreover, 
Guinotte et al. [6] predict that an increase in ocean acidification 
will reduce the ASH worldwide, decreasing the suitability of the 
ocean environment for deep-water corals. Efforts on regional and 
global scales are necessary to support the sustainable persistence of 
this fauna. 

Recent studies have drawn attention to include deep-water 
environments in marine protected areas in the Caribbean [26,44]. 
According to the present results, the assemblage of deep-water 
corals from the Caribbean are not uniform either in depth range 
or at the ecoregional scale, and the overall diversity depends 
largely on beta diversity. Therefore, an effective strategy for the 
conservation of diversity must include a broad initiative covering 
ecoregions and depth ranges. As consider the design of marine 
protected areas for deep waters in various ecoregions and define 
basin-scale standards for the exploitation of deep habitats. In 
addition, diversity may have patterns of variation at various scales 
within ecoregions, as shown by the spatial patterns in Colombia 
[24]. Intra-regional patterns, areal estimation [15,82] and 
modeling of the environmental factors that ensure the sustainabil- 
ity of coral populations at regional scale are necessary for the 
proper understanding and management of this ecosystem. 
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